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The  paper  offers  a  novel  approach  to  parameter  estimation  of  a  single-diode  solar  cell/panel  equivalent 
circuit,  based  on  analysis  of  either  technical  characteristics  supplied  by  the  manufacturer  or  user- 
obtained  experimental  I-V  curve.  The  derived  model  allows  predicting  the  solar  cell/panel  output  for 
arbitrary  environmental  conditions.  The  method  combines  a  solution  of  an  algebraic  equations  system 
with  an  optimization  algorithm.  The  main  advantage  of  the  proposed  method  is  that  in  order  to  obtain 
the  required  parameters  of  the  model,  minimal  set  of  experimental  data  is  required.  In  order  to  validate 
the  feasibility  of  the  proposed  approach,  several  solar  panels  of  different  types  from  several  manufac¬ 
turers  were  analyzed.  The  results  demonstrate  0.1 -0.5%  estimation  precision  when  the  Standard 
Operation  Conditions  data  from  the  manufacturers'  datasheets  are  employed. 
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1.  Introduction 

Modeling  of  photovoltaic  devices  and  simulating  their  behavior 
represent  significant  portion  of  today's  research  in  the  field  of 
solar  energy.  Modeling  using  an  equivalent  circuit  (lumped  para¬ 
meters  model)  is  of  particular  interest,  since  it  allows  performing 
joint  simulation  of  a  photoelectric  device  with  power  electronics 
interfaces  [1-5].  Such  simulations  allow  optimizing  the  design  of 
solar  arrays  and  power  systems  and  provide  valuable  assistance 
during  Maximum  Power  Point  Tracking  algorithms  development 
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without  costly  experiments  but  using  computer  models  only  [6-8]. 
Hence,  proper  selection  of  an  equivalent  circuit  topology  along 
with  suitable  methodology  of  extracting  optimal  set  of  equivalent 
circuit  parameters,  characterizing  a  specific  device,  is  a  very  impor¬ 
tant  part  of  modeling  process. 

Dozens  of  scientific  and  engineering  publications  covering  a 
variety  of  methods  for  extracting  the  parameters  of  equivalent 
circuit  models  were  revealed  recently.  Most  researchers  propose  a 
topology,  including  both  linear  and  nonlinear  components  [9  .  The 
linear  part  consists  of  resistive  elements  and  irradiation  dependent 
current  source.  Typically,  two  equivalent  resistors  are  included: 
series  and  shunt/parallel  resistances  [10-12].  In  some  cases,  the 
shunt  resistance  is  neglected  since  it  is  much  larger  than  the  series 
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one  [13,14].  One  or  more  ideal  diodes  typically  represent  the 
nonlinear  part  of  the  equivalent  circuit  [10,12,15,16  . 

Two  fundamental  approaches  are  used  for  extracting  the  model 
parameters  of  a  solar  cell/panel.  The  first  one  involves  multiple 
optimization  techniques,  the  most  prominent  example  of  which  is 
the  least  mean  squares  approach.  There,  the  desired  parameters 
are  obtained  by  minimizing  the  difference  between  measured 
data  and  the  modeled  I-V  curve.  The  method  requires  fluency 
in  sophisticated  algorithms  [  17—20,22—24]  and  returns  accurate 
values  of  the  desired  parameters  in  case  a  sufficiently  large 
number  of  reference  points  of  the  measured  I-V  curve  is  available 
[21,25-27]. 

Solving  systems  of  algebraic  (usually  transcendental)  equations 
is  the  second  fundamental  approach  to  parameter  extraction  of  a 
solar  cell/panel  equivalent  circuit  [28,29  .  The  amount  of  experi¬ 
mental  data  used  in  these  methods  must  match  the  number  of 
unknown  parameters.  The  equations  are  established  taking  into 
account  the  constraint  that  the  modeled  I-V  curve  should  pass 
through  three  basic  points  of  the  measured  I-V  curve.  Additional 
information  such  as  the  derivatives  of  I-V  curve  at  the  basic  points 
may  be  used  as  well  [30-34].  In  [13],  parameter  extraction  of  a 
single-diode  model  was  performed  by  solving  a  system  of  four 
algebraic  equations,  three  of  which  represent  the  I-V  curve 
passing  through  a  set  of  the  following  basic  points:  short  circuit 
current,  open  circuit  voltage  and  maximum  power.  The  fourth 
equation  used  the  fact  that  the  power  curve  derivative  is  zero  at 
the  maximum  power  point. 

The  main  advantage  of  analytical  methods  is  that  in  order  to 
obtain  the  required  parameters  of  the  equivalent  circuit,  minimal 
set  of  experimental  data  is  required  [35].  In  many  cases,  the 
required  parameters  could  be  extracted  directly  from  the  manu¬ 
facturer's  data.  The  method  in  [13]  returns  good  correspondence 
between  theoretical  and  empirical  data  for  single  crystalline 
photovoltaic  cells  and  panels  for  relatively  high  irradiation  levels. 
However,  neglecting  the  shunt  resistance  causes  significant  diver¬ 
gence  between  the  theoretical  and  real  I-V  curve  for  low  irradia¬ 
tion  levels.  Moreover,  this  assumption  is  invalid  for  modeling 
polycrystalline,  amorphous  silicon,  and  organic  solar  cells  under 
any  irradiation  level. 

In  the  present  study,  a  novel  approach  is  proposed  for  extracting 
the  single-diode  equivalent  circuit  parameters,  combining  both 
above  mentioned  fundamental  approaches,  i.e.  minimization  of 
divergences  between  the  modeled  I-V  curve  and  experimental  data 
by  fitting  and  solution  of  algebraic  equations  (partially  transcen¬ 
dental)  system.  The  approach  is  not  new  and  has  been  recently 
discussed  in  [51,52];  nevertheless,  unlike  in  the  mentioned  refer¬ 
ences,  here  the  diode  ideality  is  used  as  the  fitting  parameter  since 
its  value  resides  in  a  relatively  narrow,  well-defined  region.  In 
addition,  no  simplifying  assumptions  are  made  during  the  analysis. 
The  proposed  method  allows  reducing  the  number  of  simulta¬ 
neously  solved  equations  to  two.  The  other  two  parameters  can  be 
then  obtained  by  direct  calculation,  while  the  fifth  parameter  is 
derived  through  optimization/fitting. 

Equivalent  circuit  of  the  single  diode  model  adopted  in  the 
current  study  is  shown  in  Fig.  1.  The  model  includes  a  parallel 
combination  of  a  photogenerated  current  source  Iph,  a  diode  D, 
described  by  the  well-known  single-exponential  Shockley  equa¬ 
tion  [36  ,  a  shunt  resistance  Rsh  and  a  series  resistance  Rs.  An 
arbitrary  load  is  connected  to  the  ports  of  the  model  and  is 
characterized  by  the  current  and  voltage  IL  and  VL,  respectively. 

Note  that  an  analytical  solution  of  system  in  Fig.  1  does  not  exist  in 
general  case.  Utilizing  the  Lambert-W  function  simplifies  the  numer¬ 
ical  solution  [14,37-39].  This  function  is  the  inverse  of  the  function 
f(x)  =  xex,  i.e.  W(x)  where  W(x)  is  the  Lambert-W  func¬ 

tion.  The  solution  process  is  consequently  simplified  since  the 
Lambert-W  function  is  a  built-in  function  in  most  of  the  existing 
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Fig.  1.  The  “single-diode”  equivalent  circuit  of  the  photovoltaic  cell  with  a  load. 

engineering  software  packages,  e.g.  lambert  w  in  Matlab  and  Product- 
Log  in  Mathematica. 

A  significant  contribution  of  the  proposed  method  is  the 
possibility  of  obtaining  estimated  values  of  band  gap  energy  of 
the  cell  and  the  temperature  coefficient  of  the  photo-generating 
current.  The  value  of  the  band  gap  is  usually  assumed  constant; 
nevertheless,  its  value  varies  from  one  material  to  another  and 
strongly  depends  on  temperature  [40  .  Moreover,  since  the  single¬ 
diode  model  only  approximates  the  underlying  physical  processes, 
some  unmodeled  behavior  exists  in  any  case  and  is  reflected 
by  variations  in  band  gap  value.  The  information  regarding  the 
temperature  coefficient  of  photo-generating  current  absent  from 
datasheets  (the  temperature  coefficient  of  the  short-circuit 
current  is  usually  given).  Most  of  the  researchers  utilize  the  given 
coefficient  for  estimating  the  photo-generating  current  for  differ¬ 
ent  ambient  conditions  by  assuming  that  short-circuit  current  and 
photo-generating  current  are  equal.  Here  such  assumption  is  not 
used.  Instead,  the  temperature  coefficient  of  photo-generating 
current  is  analytically  calculated.  Equations  allowing  estimating 
the  band  gap  value  and  the  temperature  coefficient  of  photo¬ 
generating  current  require  information  regarding  temperature 
coefficients  for  both  open  circuit  voltage  and  short  circuit  current, 
which  are  possible  to  obtain  either  by  test  measurements  or  from 
manufacturer  specifications,  as  mentioned.  Thus,  only  the  elec¬ 
tron's  charge  and  the  Boltzmann  constant  are  assumed  constant  in 
the  proposed  method. 


2.  Materials  and  methods 

The  proposed  methodology  consists  of  four  main  steps.  The  first 
step  involves  obtaining  the  I-V  curve  expressions  based  on  the 
single-diode  equivalent  circuit,  also  referred  to  as  “five  parameter 
model”.  In  the  second  step,  substitution  of  the  standard  temperature 
conditions  (STC)  data  into  expressions  derived  during  the  first  step 
creates  a  set  of  four  algebraic  equations,  still  insufficient  for 
determining  the  five  parameters  of  the  chosen  equivalent  circuit. 
Thus,  the  fifth  parameter  is  determined  by  minimization  of  the 
difference  (i.e.  by  fitting)  between  the  corresponding  sets  of  mod¬ 
eled  and  experimental  I-V  curves  in  the  third  step.  Here,  the  diode 
ideality  factor  is  chosen  as  the  fitting  parameter  [37].  The  fourth  step 
involves  a  solution  of  two  additional  equations,  in  which  thermal 
coefficients  of  open  circuit  voltage  and  short  circuit  current  are 
utilized  to  determine  the  band  gap  energy  and  temperature  coeffi¬ 
cient  of  the  photo-generating  current. 

It  is  important  to  note  that  all  the  equations  are  derived  for  a 
single  while  the  manufacturer's  data  is  usually  given  for  a  panel. 
Since  a  typical  solar  panel  consists  of  multiple  serially  connected 
cells,  the  panel  voltage  is  scaled  down  by  number  of  cells  to  obtain 
single  cell  data  [41  .  Nevertheless,  in  case  a  solar  panel  consists  of 
parallel  cell  connection,  the  panel  current  should  be  similarly 
scaled. 
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3.  Theory/calculation 


3.1.  Deriving  the  set  of  equivalent  circuit  equations 


A  set  of  algebraic  equations  is  derived  according  to  Kirchhoff  s 
laws,  taking  into  consideration  the  characteristic  of  all  circuit 
elements.  The  diode  current  may  be  described  by  Shockley's 
expression  [36]  as 

/d  = /o(e(WT)  - 1),  (1) 


where  VT  [V]  is  thermal  voltage,  I0  [A]  is  reverse  bias  saturation 
current  and  rj  is  diode  ideality  factor.  The  thermal  voltage  is 
expressed  as 


Vt  =  1T  =  vt{t*)’  (2) 

where  kb= 1.380  x  10-23  [JK-1]  is  Boltzmann  constant,  q  =  1.602 
x  10“  9  [C]  is  elementary  charge  and  T  [K]  is  cell/panel  tempera¬ 
ture.  The  asterisk  (*)  here  and  thereafter  denotes  a  value  corre¬ 
sponding  to  STC,  i.e.  irradiance  of  1  kW  m~2,  AM  1.5  and  cell/panel 
temperature  of  25  °C. 

The  reverse  bias  saturation  current  depends  on  temperature  as 

[42] 


* 

o 


f  I_\3  e(Eg/qtlV*)(T-r/T) 

\T V 


where  Eg  [eV]  is  the  band  gap  energy.  According  to  [43-45],  and 
[13],  the  photocurrent  is  given  by 

Iph  =  kph^I*pll(-L+e(T-T*))  (4) 

where  e  [ K1  ]  is  the  temperature  coefficient  of  the  photocurrent, 
*F  [W  m“ 2]  is  the  irradiance  and  kph  is  the  spectral  factor  (kph  =  1 
for  AM  1.5)  [13]. 

Applying  Kirchhoff  s  laws  to  Fig.  1,  the  set  of  equivalent  circuit 
equations  consists  of  the  three  following  relations, 

Iph  ~h  —  hh  -h  =  0,  (5) 


Vd-Vl-IlRs  =  0, 


extracted  from  the  datasheet  I-V  curve  with  relatively  high 
precision  as  well. 

Substitution  of  IL  =  0  and  VL  =  V*c  into  (8)  leads  to  the  follow¬ 
ing  open  circuit  STC  relation: 


V*0c  =  (I*o+i;h)-Rsh-ri-V*T-W 


T% 


r* 


J* 
1 0 


* 


Jl  ■  vT 


Rsh.e«W^-v*A 


(10) 


Furthermore,  substitution  of  VL  =  0  and  IL  =  I*c  into  (9)  leads  to  the 
subsequent  short  circuit  STC  expression, 


* 


/*  = 


RsHRsh  n*  ,  r*  \ 


lsc 


R< 


(«h) 


R, 


W 


j* 

qV* 


RsWRsh 


a*  +  /*.)' 
Rsh\\Rs-^- 

e  1 t 


(11) 


The  (V*,/*)  voltage-current  pair  determines  the  location  of 
maximum  power  point  on  the  STC  I-V  curve, 

<  =  (12) 


Hence,  combining  (8)  and  (12),  maximum  power  point  STC 
relation  is  given  by 


K  =  do +4)  •  Rsh-%  ■  (Rs+Rsh)->i  •  v*T  •  w\ 


j* 
1 0 


•  v? 


Rsh  ■  e{ 


0 U*o-I*m+I*ph)-Rsh/r,V *) 


(13) 


The  following  condition  determines  the  extremum  of  the  power 
curve  at  h  =  Im: 


d(IL  ■  VL) 
dlL 


i,  =  n 


(14) 


Combining  of  (8)  and  (14)  leads  to 


(I%  +  I*)-Rsh-2I 


% 

m 


(Rs  +Rsh)++W  ~  c  +  'pi’ )Rsh/"  'v'} ) 


ItR 


lmJ 


sh 


r/-Vj  |  =  0  •  (15) 


1  +  W((f0  ■  Rsh/ti  ■  V?)  •  eCW-C  +';h W*/i  ^)) 


To  summarize,  a  system  of  four  transcendental  Eqs.  (10),  (11),  (13), 
and  (15)  is  now  derived  for  STC.  Note  that  the  obtained  set  of 
equations  is  still  insufficient  for  complete  solution  since  there  are 
five  unknowns. 


Vd 

Rsh 


Note  that  resistors  Rs  and  Rsh  as  well  as  the  diode  ideality  factor 
rj  are  further  assumed  to  be  temperature  independent  [13]. 
Combining  (1)  and  (5)-(7)  and  applying  the  Lambert-W  function, 
the  expression  of  I-V  curve  may  be  presented  in  the  two  following 
forms, 


V i  —  (Io  +  Iph)  '  Rsh  —  II  '  (Rs+Rsh)  —  f 1  'Vj-W 


/o 


y)  ■  Vj 


•  e^°  -E+IphV  Rsh/bl  '  VtJ) 


and 


j  do  +  Iph) '  Rsh  —  Vl  ryVr 


Io 


R,R 


sh 


Rs  +  Rsh 


Rs 


rj-VT  ( Rs  +  Rsh ) 


,(Rsh  ■  ((Io  +  Iph)  ‘  Us  +  VL)/r]  ■  VT  ■  (Rs  +Rsfi))^ 

(9) 


3.2.  Substitution  of  STC  data  into  the  I-V  curve  expressions 

Expressions  (8)  and  (9)  allow  to  carry  out  the  calculations  for 
determination  of  all  the  five  unknown  equivalent  circuit  para¬ 
meters,  namely  Iph,I0,Rs,Rsh,  and?/.  The  STC  data  is  typically 
provided  by  the  manufacturer  and  is  utilized  as  follows.  The  short 
circuit  current  /*c,  the  open  circuit  voltage  V*c  and  the  maximum 
power  point  voltage-current  pair  (V*  ,/* )  of  the  STC  I-V  curve  are 
referred  to  as  basic  points.  The  manufacturer  usually  supplies 
these  data  in  a  tabular  form;  moreover,  these  points  can  be  easily 


3.3.  Solving  the  system  of  equations 


A  possible  way  to  solve  the  derived  system  of  equations  is 
treating  one  of  the  unknowns  as  an  independent  parameter,  which 
may  be  determined  by  fulfillment  of  an  additional  condition, 
namely  minimization  of  deviation  between  the  theoretical  and 
measured  I-V  curves.  Each  one  of  the  five  unknowns  may  be 
selected  as  an  independent  one,  though  the  ideality  factor  is 
probably  the  most  convenient  since  its  value  resides  in  a  known 
bounded  range  for  any  type  of  photovoltaic  cells  (2  <rj<  2  for 
silicon  cells/panels)  [37]. 

Combining  (1),  (5)  and  (6)  for  open  and  short  circuit  STC 
conditions  results  in 


!ph  •  Rsh  -  Kc 

0  Rsh  ■  (e+'oc/”  ■  v’r)  -  i) 

And 


(16) 


Rsh  +  Rs  4(Rsft+Rs)-C  _  i 

Ph  SC  Rsh  Rsh  e(Voc/r,V*)_e(l*cRs/„VTy  y  > 

respectively.  Substitution  of  (16)  and  (17)  into  (13)  and  (15)  results 
in  the  system  of  two  transcendental  equations  with  three 
unknowns  Rs ,  Rsh  and  rj,  where  rj  is  an  independent  parameter, 
residing  within  a  known  range.  The  system  of  these  two  trans¬ 
cendental  equations  is  solvable  only  by  numerical  methods.  Each 
value  of ;/  creates  a  specific  set  of  ( Rs ,  Rsh),  for  which  /J  and  Iph  may 
be  found  from  (16)  and  (17).  Each  full  set  of  five  parameters 
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determines  a  specific  l-V  curve  that  passes  through  the  three  basic 
points  and  possesses  a  maximum  power  point  at  (Jm,  Vm)  as  shown 
in  Fig.  2.  Certain  deviation  exists  between  the  curves  except  at  the 
three  basic  points,  where  the  curves  must  coincide.  The  correct 
value  of  rj  corresponds  to  the  modeled  curve  with  minimum 
deviation  from  the  datasheet/measured  l-V  curve  at  STC. 

The  minimization  may  be  performed  by  visual  comparison 
between  modeled  curves  and  the  measured  one  as  well  as  by  the 
simplest  error  minimization  algorithm.  When  a  measured  l-V 
curve  is  unavailable  in  raw  data  form  and  a  datasheet  picture  only 
is  available  (which  is  a  typical  case),  one  can  use  digitizer  software, 
e.g.  [45],  which  helps  digitizing  scanned  plots  into  a  raw  data  with 
fair  precision. 


3.4.  Additional  set  of  equations  for  thermal  parameters  calculation 

The  first  three  steps  allow  a  fairly  accurate  reproduction  of  the 
STC  l-V  curve.  The  final  step  is  estimating  l-V  curves  for  arbitrary 
environmental  conditions  using  temperature  coefficients.  The 
expression  for  the  open  circuit  voltage  temperature  coefficient 
(TCVoc)  can  be  obtained  by  combining  the  derivative  of  (8)  with 
(2)-(4)  and  substituting  T=T , 


where 


(23) 


The  values  of  Eg  and  e  can  be  then  obtained  by  solving  (20) 
and  (22). 


4.  Results  and  discussion 

The  proposed  method  was  first  applied  to  Isofoton  IS  -  160 
mono-crystalline  solar  panel  46].  The  seven  parameters  of  the 
model  ( Rs ,  Rsh,  Iq,  I*h,  rj,  Eg  and  e)  were  obtained  in  accordance 
with  the  four  steps  of  the  presented  approach.  Theoretical  l-V 
curves  derived  using  the  proposed  method  for  different  environ¬ 
mental  conditions  were  compared  to  experimental  l-V  data 
provided  by  manufacturer.  All  the  equivalent  circuit  parameters 
were  obtained  from  the  manufacturer's  data  only.  The  results  were 
obtained  for  the  same  environmental  conditions  as  manufacturers' 
and  transferred  into  transparent  l-V  curve  plots.  These  plots  were 
superimposed  with  the  manufacturers'  datasheets  by  correspond¬ 
ing  scaling  along  each  axis. 

The  specifications  of  the  IS  -  160  panel  and  the  manufacturer's 
l-V  curve  are  shown  in  Table  1  and  Fig.  3,  respectively.  Note  that 


dVL 

dT 


T  =  T* 


EgI*0Rsh  +  qRsh( 3/J  +  i;h  T*£)W*t  -  n V*A(Eg  -  qRsh(I*0  +  I*ph  -  IL) +3 qqV*  +  qr,V*A) 

r,V*(A  +  \)qT* 


(18) 


where 


A  =  w  (Il^lle(RsH(l,o+i;h-k)/nv,T)) 

\Wt  ) 


(19) 


The  TCVoc  is  obtained  from  (18)  by  substituting  IL= 0  as 

EgI*0Rsh  +  qRsh(3I*0  +I*phT*e)r,V*  -  rjV*B(Eg  -  qRsh(I*0 +  l*ph)  +  3qr,V*  +  qrjV*B) 


TCVoc  = 


rjqT*V*(B+ 1) 


(20) 


where 


B  =  w(  !°^!ie(Rih(1o +  ) . 


* 


Y]V  t 


(21) 


Similarly,  the  expression  for  the  short  circuit  current  temperature 
coefficient  (TC/sc)  may  be  derived  from  combining  the  derivative  of 
(9)  with  (2)-(4)  and  substituting  T=T* , 


the  data,  provided  by  Isofoton  in  tabular  form  and  the  datasheet 
l-V  curve  slightly  differ  (common  to  different  solar  panels  man¬ 
ufacturers).  The  deviation  remains  within  the  tolerance  ( +  5%) 
declared  by  the  manufacturer. 

Since  the  proposed  method  is  very  sensitive  to  input  data, 
graphical  data  was  used  here  rather  than  tabular  one.  The  estima¬ 
tion  procedure  was  carried  out  using  Mathematica  software  [47]. 
The  estimation  results  for  STC  are  summarized  in  Table  2. 

The  estimated  parameters,  given  in  Table  2,  were  used  for 
creating  modeled  l-V  and  P-V  curves  for  STC  (Fig.  4)  and  different 
temperatures  (Fig.  5).  Good  agreement  may  be  observed,  char¬ 
acterized  by  deviations  of  0.1 -0.5%. 

Furthermore,  the  proposed  method  was  applied  to  three 
additional  solar  panels,  specified  in  Tables  3-5.  The  estimation 
results  are  summarized  in  Tables  6-8.  The  corresponding  l-V  and 


Trj  _  dIL 
1  -  dT 


T  =  T* 


RsllRsh 

r,V*qRsT*(C+ 1) 


( V  o  +>/V?<J(3/S  +I*phT*e)  +  q(I*0  +  i;h) 


Eg-3gV*Tq-qV^qC  qV*  Q \ 

RshWRs  RsRsh  ) 


(22) 


Fig.  2.  The  family  of  l-V  curves  with  different  ideality  factors  passing  through  the 
basic  points 


Table  1 

Specification  of  the  IS  -  160  solar  panel  [46]. 


Electrical  characteristics,  STC 

From  datasheet 
table 

From  datasheet 

curve 

Number  of  series-connected 

72 

72 

cells  ns 

Maximum  power  Pm,  [W] 

160  (  ± 3%) 

155.67 

Maximum  power  voltage  Vm, 

35 

35.2 

[V] 

Maximum  power  current  Im,  [A] 

4.57 

4.42 

Shot  circuit  current  Isc,  [A] 

4.9 

4.81 

Open  circuit  voltage  Voc,  [V] 

43.8 

43.2 

Temperature  coefficient  of  Voc, 

-0.378 

-0.353 

%/K 

Temperature  coefficient  of  Isc, 

0.0254 

0.0254 

%/K 
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Fig.  3.  The  STC  I-V  curve  of  IS  -  160  panel,  digitized  from  the  datasheet. 


Table  2 

The  estimated  STC  parameters  for  a  single  cell  of  the  IS-160  panel. 


Parameter 

Estimated  value 

Ideality  factor,  rj 

1.0 

Series  resistance  Rs,  [mQ] 

8.0 

Shunt  resistance  Rsh,  [X2] 

3.25 

Photocurrent  Iph,  [A] 

4.824 

Reverse  bias  saturation  current  /0,  [pA] 

0.33 

Band  gap  energy  Eg,  [eV] 

1.16 

Temperature  coefficient  of  Iph,  [pA/K] 

254 
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Fig.  4.  The  I-V  and  power  curves  of  IS-160  at  STC.  Simulation  results  (dotted)  are 
added  to  the  manufacturer's  data  (solid). 


/-  V  Curve  variation  according  to  temperature 
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Fig.  5.  The  I-V  and  power  curves  of  IS-160  for  different  temperatures.  Simulation 
results  (dotted)  are  added  to  the  manufacturer's  data  (solid). 


Table  3 

Specification  of  the  Suntech  CNPV-280P  solar  panel  [48]. 


Electrical  characteristics,  STC 

From  datasheet  I-V  curve 

Number  of  series-connected  cells  ns 

72 

maximum  power  Pm,  [W] 

280 

Maximum  power  voltage  Vm,  [V] 

36.9 

Maximum  power  current  Im,  [A] 

7.6 

Shot  circuit  current  Isc,  [A] 

8.2 

Open  circuit  voltage  VOCl  [V] 

44.6 

Temperature  coefficient  of  Voc,  %/I< 

-0.3 

Temperature  coefficient  of  Isc,  %/K 

0.05 

Table  4 

Specification  of  the  Schuco  SF-160-24-M175  solar  panel  [49]. 

Electrical  characteristics,  STC 

From  datasheet  I-V  curve 

Number  of  series-connected  cells  ns 

72 

Maximum  power  Pm,  [W] 

175 

Maximum  power  voltage  Vm,  [V] 

36 

Maximum  power  current  Im,  [A] 

4.86 

Shot  circuit  current  Isc,  [A] 

5.2 

Open  circuit  voltage  Voc,  [V] 

44.8 

Temperature  coefficient  of  Voc,  %/K 

-0.48 

Temperature  coefficient  of  Isc,  %/K 

0.04 

Table  5 

Specification  of  the  Kyocera  KC200GT  solar  panel  [50]. 

Electrical  characteristics,  STC 

From  datasheet  I-V  curve 

Number  of  series-connected  cells  ns 

54 

Maximum  power  Pm,  [W] 

200 

Maximum  power  voltage  Vm,  [V] 

26.3 

Maximum  power  current  Im,  [A] 

7.61 

Shot  circuit  current  Isc,  [A] 

8.21 

Open  circuit  voltage  Voc,  [V] 

32.9 

Temperature  coefficient  of  Voc,  %/K 

-0.374 

Temperature  coefficient  of  Isc,  %/K 

0.039 

Table  6 

The  estimated  STC  parameters  for  a  single  cell  of  the  CNPV-280P  panel. 

Parameter 

Estimated  value 

Ideality  factor,  rj 

1.0 

Series  resistance  Rs,  [mO] 

4.2 

Shunt  resistance  Rsh,  [Q] 

4.3 

Photocurrent  Iph,  [A] 

8.19 

Reverse  bias  saturation  current  /0,  [nA] 

0.27 

Band  gap  energy  Eg,  [eV] 

1.19 

Temperature  coefficient  of  Iph,  [pA/K] 

420 

power  curves  of  these  panels  along  with  the  manufacturers' 
supplied  data  are  shown  in  Figs.  6-8.  It  can  be  concluded  that  the 
estimated  curves  nearly  coincide  with  the  manufacturer  provided 
curves,  indicating  the  feasibility  and  high  precision  of  the  proposed 
method. 


5.  Conclusions 

An  improved  method  for  determining  the  parameters  of  single¬ 
diode  photovoltaic  cell  equivalent  circuit  was  developed  in  the 
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Table  7 

The  estimated  STC  parameters  for  a  single  cell  of  the  SF-160-24-M175  panel. 


Parameter 

Estimated  value 

Ideality  factor,  rj 

1.2 

Series  resistance  Rs,  [mCi] 

7.1 

Shunt  resistance  Rsh,  [Q] 

10.2 

Photocurrent  Iph,  [A] 

5.18 

Reverse  bias  saturation  current  /0,  [nA] 

10.2 

Band  gap  energy  Eg,  [eV] 

1.16 

Temperature  coefficient  of  Iph,  [pA/I<] 

300 

Table  8 

The  estimated  STC  parameters  for  a  single  cell  of  the  KC200GT  panel. 

Parameter 

Estimated  value 

Ideality  factor,  rj 

1.05 

Series  resistance  Rs,  [m^] 

4.85 

Shunt  resistance  Rsh,  [Q] 

3.6 

Photocurrent  Iph,  [A] 

8.2 

Reverse  bias  saturation  current  /0,  [nA] 

1.05 

Band  gap  energy  Eg,  [eV] 

1.22 

Temperature  coefficient  of  Iph,  [pA/K] 

388 

presented  work.  The  set  of  experimental  data  required  for  the 
extraction  of  model  parameters  included  manufacturer  provided 
STC  data  only. 

In  order  to  calculate  the  five  unknown  parameters  for  STC 
conditions,  the  method  combined  a  solution  of  algebraic  equations 
system  to  obtain  four  unknown  parameters;  whereas  the  fifth 
parameter  was  obtained  by  minimization  of  the  divergence 
between  the  modeled  and  experimental  I-V  curves.  Two  addi¬ 
tional  “thermal”  parameters  were  derived  from  the  solution  of 
another  system  of  algebraic  equations  based  on  the  solar  cell/panel 
temperature  coefficients.  A  number  of  I-V  curve  points  besides  the 
three  basic  points  (short  and  open  circuits,  maximum  power)  were 
required  to  obtain  an  accurate  estimation  of  the  ideality  factor. 
One  can  obtain  these  points  directly  from  the  STC  I-V  curve 
available  in  the  cell/panel  datasheet.  For  this  purpose,  dedicated 
digitizing  software  was  applied. 

Since  the  accuracy  of  the  model  parameters  strongly  depends 
on  the  accuracy  of  the  input  data,  the  accuracy  of  the  data  offered 
by  the  manufacturer  (typically  3%)  is  not  always  enough.  The 
precision  has  to  be  in  the  1-1.5%  range  for  obtaining  correct  model 
parameters.  High  accuracy  of  the  temperature  coefficients  is  also 
necessary.  Correction  of  manufacturer  provided  is  often  necessary 
and  may  be  performed  by  deriving  these  parameters  from  the 
manufacturer  provided  graphical  data. 
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Fig.  6.  The  I-V  curves  of  CNPV-280P  for  different  irradiance  levels  and  temperatures.  Simulation  results  (dotted)  are  added  to  the  manufacturer's  data  (solid). 


Fig.  7.  The  I-V  curves  of  SF  160-24-M  for  different  irradiance  levels  and  temperatures.  Simulation  results  (dotted)  are  added  to  the  manufacturer's  data  (solid). 
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Fig.  8.  The  I-V  curves  of  KC200GT  for  different  irradiance  levels  and  temperatures.  Simulation  results  (dotted)  are  added  to  the  manufacturer's  data  (solid). 


The  proposed  method  was  verified  by  comparison  between 
modeled  and  datasheet  I-V  curves  of  several  off  the  shelf  panels  of 
leading  manufacturers.  The  results  demonstrated  average  devia¬ 
tion  of  0.1 -0.5%,  indicating  the  proposed  method  feasibility  and 
high  accuracy. 
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